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Abstract. 

Photoassociation (PA) spectrum in the presence of a magnetic Feshbach resonance 
is analyzed. Nonperturbative solution of the problem yields analytical expressions 
for PA linewidth and shift which are applicable for arbitrary PA laser intensity and 
magnetic field tuning of Feshbach Resonance. We show that by tuning magnetic field 
close to Fano minimum, it is possible to suppress power broadening at increased 
laser intensities. This occurs due to quantum interference of PA transitions from 
unperturbed and perturbed continuum. Line narrowing at high laser intensities is 
accompanied by large spectral shifts. We briefly discuss important consequences of 
line narrowing in cold collisions. 

PACS numbers: 34.50.Cx, 34.80.Dp, 32.70.Jz, 34.80.Pa 
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1. Introduction 

Over the years, photoassociation and Feshbach resonance have become important tools 
in manipulation of ultracold collisions. Recent experimental [H [21 E] and theoretical 
[H El O [7] works on photoassociation (PA) near a magnetic field Feshbach resonance 
(MFR) [S] give rise to the exciting possibilities of coherent control of atom-molecule 
conversion. In a recent experiment, Junker et al. [1] have demonstrated asymmetric 
spectral line shape and saturation in PA in the presence of MFR. Asymmetric line 
shape is a characteristic feature of Fano effect [9J which arises in different areas such as 
atomic [lOj, particle [11] and condensed matter physics [12j. When a pair of colliding 
atoms under the influence of Feshbach resonance are photoassociated into an excited 
molecular state, PA transitions can occur in two competing pathways. The presence 
of a Feshbach resonance largely perturbs the continuum states of colliding atoms. As 
a result, continuum states get hybridised with one or more bound states embedded in 
the continuum. PA transitions from unperturbed and perturbed continuum states can 
interfere resulting in asymmetric Fano line shape. The unique aspect of MFR induced 
Fano effect in PA is that the continuum-bound coherence can be controlled by tuning 
the magnetic fleld. Feshbach resonances and photoassociation have a common feature: 
Both are the effects of continuum-bound interacting systems. In the case of Feshbach 
resonance, an initial continuum state is coupled to a bound state embedded in the 
continuum by hyperflne interaction of two atoms. In the case of PA, the continuum 
state gets coupled to an excited diatomic molecular state by a single photon. Both can 
be treated within the framework of Fano's theory. In fact, based on Fano's method, 
theory of one- and two-photon PA have been developed by Bohn and Julienne p!3l [H] 
Two-photon PA [15] has been shown to lead to quantum interference. Furthermore, 
quantum interference has been demonstrated in coherent formation of molecules jl6j. 

Apart from MFR-induced modiflcation, PA under intense laser flelds (strong- 
coupling regime) can further modify the continuum states due to strong continuum- 
bound dipole couphng. Fano effect in the strong fleld regime [171 [181 HI] leads to 
"confluence of bound-free coherences" [T7] which can result in a number of effects 
such as line narrowing [20], Autler-Townes splitting [17] and nonlinear Fano effect |21j . 
Continuum-bound coherences in strongly interacting ultracold atoms can lead a number 
of profound physical effects which may not be observed in the weak coupling regime. A 
decade ago, Vuletic et al. [22] experimentally observed line narrowing in the spectrum 
of trap loss of atoms due to a tunable Feshbach resonance. In Fano's Theory [9], a 
continuum interacting with a bound state is exactly diagonalized leading to a "dressed" 
continuum state. In linear Fano Theory, optical transition matrix element between 
this dressed continuum and any other bound state is obtained perturbatively by the 
use of Fermi Golden rule. In the strong coupling regime, optical coupling of the dressed 
continuum with any other state needs to be treated nonperturbatively either by "double 
diagonalization" technique [T^ or by other diagonalization techniques such as used in 
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Figure 1. A schematic diagram showing the coupling between the bound state 
(magenta) | 2) and and continuum (green) | E) with the excited bound state (blue) 
I 3) via same laser Lpa{ magenta and green double-arrow vertical lines). Red dashed 
line indicates the hyperfine coupling between the open and closed channel. 



Here we explore the possibility of suppression of power-broadening in strong- 
coupling PA by manipulation of continuum-bound coherences with Feshbach resonance. 
As the simplest possible model, we consider optically coupled two bound states 
interacting with a common continuum as shown in figure 1. We demonstrate that by 
tuning the magnetic field close to Fano minimum where excitation probability vanishes 
[9], it is possible to obtain line narrowing in PA spectrum with large shifts at high laser 
intensities. Large light-shift with narrow linewidth may be useful in efficient tuning of 
elastic scattering length by optical means [211 [2S1 126] . 

The paper is organised in the following way. In the following section, we discuss in 
brief the formulation of the problem. In section 3, we present the analytical results on 
line narrowing and enhencement of shifts. In section 4, we discuss the numerical results 
and the paper is concluded in section 5. 

2. Formulation of the problem 

The dressed state of a system of two colliding atoms interacting with a PA laser in the 
presence of MFR can be written as 

(1) 

where r is the relative coordinate of the two atoms, E is an energy eigenvalue, 
I 1(2)) represents the internal electronic states of open(closed) channel | 1(2)) and 



r 



/ dE'bE'^E'ir) I 1) + ^ $ 



i=2,3 
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I 3) denotes the electronic state of the excited molecule. {i = 2,3) and "^e' are 
the wave functions of perturbed bound and continuum states, respectively. Here Be' 
is density of states of the unperturbed continuum. In writing the above equation, 
we have assumed that PA laser can couple only a particular ro-vibrational state of 
the excited molecule; and the bound state $2 has zero angular momentum. Let the 
hyperfine spin coupling between the channels 1 and 2 be denoted by V{r). Let Qi{r) 
and ^l2{r) represent the molecular Rabi couplings of the excited state | 3) with the 
ground states | 1) and | 2), respectively. In the absence of these three couplings, 
let the unperturbed bound states be denoted by 03 (r) and 02 (?^) with bound state 
energies £^3 and E2 respectively; and the unperturbed continuum states by ipE' with 
asymptotic collision energy E'. With the use of these unperturbed solutions, we 
construct three Green's functions GE{r,r'), G2{r,r') and (j3(r, r') which correspond 
to the channels 1, 2 and 3, respectively. The continuum Green's function GE{r,r') 
can be written as GE{r,r') = —T[ip'^^^{r^)ip~^{r^)^ where implies either r or r' 

whichever is smaller (greater) than the other. Here iPe^t) — V'e'^ + ^V'^;'^^ where ip^E^ and 
ip^E^ represent regular and irregular scattering functions, respectively. Asymptotically, 
^o,reg^^^ ~ jo COS rjo — siu i^q and i'^E^'^ir) ~ — (nocosr^o + josin?7o), where jo and no 
are the spherical Bessel and Neumann functions for partial wave £ = (s-wave) and 
?7o is the s-wave phase shift in the absence of laser and magnetic field couplings. The 
other two Green's functions correspond to bound state solutions and are of the form 
G-i{r,r') = - j^s^E_E^ Mr)Mr') and G2(r,r') = -;g^02(r)02(r'). Here 6 = ujl ~ uja 
is the laser-atom detuning, the bound state energy is measured from the excited 
state threshold Hua- In order to include the spontaneous emission of the excited state, 
following the prescription given by Bohn and Julienne [13], we introduce in our model 
an artificial channel with channel state ipart and its coupling Vart with excited state | 3); 
and thereby identify spontaneous linewidth h'j = 27i\{ipart \ Vart \ 03) 

Next, let us introduce a low energy dimensionless interaction parameter 

P{k) = - cot r]res{k) ~ (kares)''^ + ^^e^ (2) 

where rjresik) is the Feshbach resonance phase shift and re is related to Feshbach 
resonance linewidth ~ kG by r^ = 2h/{fiG) where C is a constant and fi is the 
reduced mass of the two atoms, ares is related to the applied magnetic field B by 
O'res = ~ B~Bo P^ - whcrc abg is the background scattering length. Bo is the resonance 
magnetic field, A is a parameter which depends on F^ and magnetic moments of the 
atoms at Bq. Having done so, we can obtain the solutions in the form $3(r) = ApA(f)3{r) 
and ^2{r) = Acc(p2{r). The explicit form of ApA [23] is given by 

. ^ exp{ir]o){l3 + q)7rnsE 

iP + z){Ap + 2/1(7 + rpA)/2} - hTpAiq - ^ ^ 

where Ap = h6 + E — {E3 + E^f^^j-^), E^j^^^^ is the energy shift in the absence of Feshbach 
resonance, VL^e = J dr(l)^{r)VLi{r-)ip'^^^ is the continuum-bound molecular dipole coupling 
and TpA is the stimulated linewidth of PA given by Fpa = 2'k\^I^e\'^- Here q is Fano's 
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asymmetry parameter defined by 

^32 + V32 . ■ s 

^ = T/ 

where (^32 = / (ir03(r)fi2(^)02 is the bound-bound Rabi coupling and V23 = 
/ J dr'dr(j)2{r)V{r)Re[GE{r,r')]Q{r')(f)3{r') represents an effective interaction between 
the two bound states mediated through their couphngs with the common continuum. 
Note that q is independent of laser intensity. In the limit k ^ 0, both Q32 and V32 become 
energy-independent while fl3EV2E ~ k. Thus at low energy, g ~ 1/k. Detailed derivation 
of the dressed state of is given elsewhere [23], but for completeness we reproduce the 
derivation in Appendix A. 

3. Analytical results 

The loss of atoms due to decay of excited state into the decay channels (which is modeled 
as an artificial scattering channel) is described by the PA rate 

KpA =< lii^art I Vart | ^) T > 

1 r ,E.^^Jt£^^^tpJIfl^^ (5) 



where (■ ■ ■) means thermal averaging over the collision energy and h6p = h6 — {E^ + 
^shift ~ ^th) is the detuning parameter. Here Eth is the threshold of the open channel. 
Here Qt = {'2'n'fiKBT/h'^), Ek = E — Eth and Kb is the Boltzmann constant. F is the 
linewidth (in the presence of MFR) given by . 

T = fiq,P)T,= ^-^^TpA. (6) 

The extra shift caused by MFR is given by 

'{q^-l)(3-2q 



AE = - 
2 



^2 + 1 ^'PA. (7) 

(EI) shows that F depends on a nonlinear function (/(g, /?)) of q and p. Note that when 
f3 ±00, that is, far away from MFR, F TpA- It is to be further noted that when 
j3 = — g, we have F = and = at which PA ceases to occur. The Fano minimum 
is given by Pmin = —q or equivalently the corresponding magnetic field Bmin- Thus 
F can be made arbitrarily small by tuning P close to —q. It is possible to suppress 
power-broadening at increased laser intensities by tuning the magnetic field B close to 

Bmin • 

Next, we discuss the weak-coupling limit of ([6]) and (JTj) when laser intensity is low. 
For this we first find the dressed continuum state in the limit TpA 0. As Fp^ ^ 0, 
ApA ^ 0, Acc can be expressed as 



^CC = -\ -TT exp i{T]o + r]res) siu r]res (8) 

y vri r 

and ipE^r) becomes 

ipE{r) = exp i{r]o + r]res) [iPe''^^ cos rj^es + i's'"'"' sin rjres] (9) 
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Figure 2. Subplots (a) and (b) show T/TpA Vs (kares) ^ and subplots (c) and (d) 
exihibit AE/TpA Vs (kares)^^- Plots are for different values q as indicated in the 
figure, k is defined as = KgT 




Figure 3. Solid lines represent the linewidth F (in unit of Fp^) given by the expression 
© as a function of {kaj.es)~^ for q = —5.0 (upper panel) and q = —0.1 (lower panel). 
Dotted lines represent the low coupling expression of pT|) for Ci — 6 and 6*2 = 1 
(upper panel); and Ci =0.1 and C2 = (lower panel). The values of Ci and C2 are 
so chosen such that in the limit Ures 0, F — > F""^"*^. 
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So in the limit of T pa ^ 0, at low energy the state | ^'_b) of ([T]) reduces to 

I ^i?)o = liAccMr) I 2) + / hEii^E'{r)dE' \ 1)]. (10) 

Taking h'^ = S{E — E'), the stimulated linewidth T"'^'*'^ in the weak coupling limit is 
given by the Fermi golden rule expression 

yweak ^ 2,r| I r0yj(r)(3 | ^]l(r) | ^Ehdrl"" 

= Tpa\1 + Ci smr]res + C2tan?7resP (11) 

where Ci = nj^/n^E and C2 = {-,/J^)n32/^3E ■ Here n^'^ = (3 | fii | t/-^^""). The 
expression ( ITTi) is in agreement with (6) of [5]. When rjres 7r/2, T'"^"'^ diverges and 
hence ( ITTj) is not valid near T^^es = vr/2. In other words, (ITTI) is not applicable close to 
Feshbach resonance. 

Here we show how the exact expressions of F and AE as given by ([6]) and ([7]), 
respectively, do enable us to realize the possibility of suppression of power-broadening 
with enhenced light shift at increased laser intensities. It is clear from ([7]) that AE 
goes to zero as /5 ±00. From ([6]) and ([7]), we notice that if laser intensity / is 
increased by a factor M/, to suppress power-broadening /3 is to be changed to /3' such 
that f{q,P') = Mf^f{q,P). For /5 ~ — g we have AE ^ —qTpA/2 which is proportional 
to /. Hence when line broadening is suppressed by the tunability of MFR, the total 
shift E shift = E^hift + remains proportional to I. The shift E^hij^^ is negative in 
the low energy regime. However, the total shift Eghift in the presence of MFR can be 
positive or negative depending on the values of q and /3. 

Finally we prove that line narrowing in one-photon PA is not possible in the absence 
of coupling between open and closed channel. It can be noticed that the PA laser can be 
tuned either near continuum-bound frequency in which case hSp ~ Ef^ or near bound- 
bound transition frequency {E3 — E2)/h in which case —Tidp = {E2 — Eth) = E^es- Here 
E^ = E^ + -^s/ii/t and E2 = E2 + AE2, AE2 being the shift of closed channel bound 
state due to its coupling V with the open channel. In the limit F^ 0, 

KpA ^ r2 (12) 

[{E-E,- EL/J + S| + '^'(^ + rPA)V4 

which is in agreement with the expression of l^i^p of [H] if we identify Aiand A2 of 
[11] with —h6 and E2, respectively. In our case there is only one laser coupling between 
the continuum and the excited bound state and also between the two bound states. It 
is clear from the above expression that in the absence of coupling between open and 
closed channels the narrowing of PA linewidth is not be possible. 

4. Numerical results and discussion 

For numerical illustration, we consider a model system of two ground-state (S1/2) Li 
atoms undergoing PA from the ground molecular configuration '^S^ to the vibrational 
state f = 83 of the excited molecular configuration 1 which correlates asymptotically 
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300 



Figure 4. KpA in cm"^ scc^^ Vs. detuning 6p in MHz. Each pair of dashed and sohd 
curves are obtained for TpA = 1.0 MHz (left pair), TpA = 10.0 MHz (middle pair) 
and TpA ~ 25.0 MHz (right pair) for the fixed q — —6.36. For solid curves, magnetic 
fields are B = 705.00 G (left), B = 708.56 G (middle) and B = 709.12 G (right). 
For dash-dotted curves, these are B = 713.19 G (left), S = 711.09 G (middle) and 
B — 710.68 G (right). The magnetic fields are so chosen such that the linewidth F 
remains fixed at 0.04 MHz. 
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Figure 5. The locus of B and TpA for which the linewidth F remains fixed at 0.04 
MHz. The lower inset shows the variation of AE (in MHz) against TpA (in MHz) and 
the upper inset exhibits the variation ares (in nm) against B (in Gauss) at the fixed 
F = 0.04 MHz. 
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to 2Si/2 + 2Pi/2 free atoms [2Sl[2n]. AH the relevant parameters 7, E^^^-j^, A, abg and 
are estimated from [28], [30] and [31]. In figure 2, we have plotted T/TpA and AE/TpA 
against 1/kares for positive and negative q values. The maximum and minimum values 
of linewidth would be observed for /? = 1/g and /? = —q, respectively. The magnitude of 
the change in shift due to PA in the presence of MFR is significant near f3 = —q. Figure 
3 clearly shows that the stimulated linewidth in the weak coupling limit, represented 
by dashed lines deviate appreciably from nonperturbative results as shown by solid 
lines. The deviations are the most prominent in the region {kares)~^ — {r]res — 7r/2). 
Furthermore, for lower q values these two results deviate most significantly. Figure 4 
illustrates how to suppress power-broadening by the appropriate tuning of magnetic 
field near Bmin and thereby to keep the total linewidth close to the natural linewidth. 
There are two values of l3{B) and correspondingly two values of ares where the linewidth 
F can be kept fixed at a small value at an increased laser intensity. In figure 5, we show 
how to vary Fp^ (or laser intensity) and the magnetic field in order to keep F fixed 
at 0.04 MHz which is much smaller than the natural linewidth 7 although TpA can be 
many orders of magnitude higher than 7. The lower inset in figure 5 shows that the 
extra shift AE can exceed 7 by many orders of magnitude while power-broadening is 
suppressed. 

5. Conclusion 

In conclusion, we have demonstrated that linewidth of photoassociation spectrum can be 
narrowed down close to the natural linewidth by making use of tunability of a magnetic 
Feshbach resonance. Experimentally, the line narrowing effect may be inferred from 
PA spectra near Bmin at a high intensity. This enhencement of the life time of excited 
molecular state may be beneficial for population transfer from ground state coUisional 
continuum to ground molecular state by two-photon Raman-type PA. Furthermore, 
narrow linewidth with large shift will be useful for efficient manipulation of scattering 
length by optical Feshbach resonance [211 [25] . This will be particularly important for 
altering scattering amplitude of higher partial waves [26] . 
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Appendix-A 

The coupled differential equations for the model are given by 



2/i dr'^ 



+ Bj{r) $3+ [Ve{r)-h6~E-zh-f/2]<l>3 



(A.l) 
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2/i dr'^ 



+ V2{r) - E 



$2 = -n;<^s - v*x, 



(A.2) 



2/i dr'^ 



+ Vi{r)-E 



(A.3) 



Using Green's functions G3{r,r'), we find $3 = Apa4>3, where 

-Mr)- 



A 



PA 



h6 + E - E3 + ih-f/2 
In a similar way, using G2{r,r') we can write $2 = Acc'4'2, where 

[ApAni, + J dr'x{r')n,{r')Ur')] 



(A.4) 



A 



cc 



E — E2 

where Q32 = J '^'"03(^)^2(^)02(^)- The continuum state can be written as x(r 
/ dE'bE'ipE'ir). Thus we obtain 

h'' d^ 



(A.5) 



2/i dr'^ 



^E'ir) + [V,ir) - E]^E'ir) = -niir)ApAMr) 



{V2E + ApAn2) 



where 



A 



E — E2 

^3EiE — E2) + f^32^2£; 



V{r)Ur) 



PA 



(A.6) 



(A.7) 



viE - E2) - \n;2\' ■ 

Here V = HS + E - E3 + ih-f/2, V2E = J dr(j)2{r)V{r)iljE'{r) and Q3E = 
J dr(l)3{r)VLi{r)ipEi{r) . Using the Green's function GE{r,r'), from equaion (]A.6P we 
obtain 



ijE' = exp(^r/o)V^^r^ + / dr'GE{r,r') \ni{r')ApAUr' 



+ 



'^2E + ^PA^32 



V(r')02(r') 



(Ai 



{E - E2) 

Now using this solution, we can calculate the probability amplitude of excitation in the 
following form 

exp(ir/o) [n^E + GV^^) 



ApA 



V 



1^32^ 
E-E2 



{B, + GBf) 



where 



and 



V2E = J rfr02(r)\/(r)V'^r^(r) 
^3E = / rfr03(r)ili(r)?/^^r^(r). 



(A.9) 

(A.IO) 
(A.ll) 
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The other parameters here are 

[E - E2 - AE2 + tTT\V2E\^ 



^ ^ [^32 + V32 - inV2E^3E] 



and 



where 



^P- ^ shift -^T^\^hE\ -^^32 - E2) lA.ldj 

Bf = V,2 - lW2E^3E - ^32^^^^;^^^^^^ (A.14) 



ELft = J J dr'drUrWiir)Re[GE{r',r)]n,{r')(Ps{r'), (A.15) 
V'32 = / J dr'dr(l)2{r)V{r)Re[GE{r',r)]ni{r')(j)3{r') (A.16) 

AE2= [ [ dr'dr(j)2{r)V*{r)Re[GE{r',r)]V{r')(j)2{r'). (A.17) 



and 



Here we introduce two new parameters (3 = {E — E2 — AE2)/{TiTr/2) and q = 
{fl32 + V32)/(7rfi3£;V2£;). Now after some algebra we can express (\A.9\i in the following 
form 

7 ^ exp{ir]o){P + q)nn3E , . 

{(3 + t){A, + zn{^ + TpA)/2}-nTpA{q-ty/2 ^ ■ ' 
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